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Abstract 

We analyze the Rice-Ruina state and rate dependent friction model. 
The system consists of one or two blocks driven by springs with constant 
velocity on a dry, rough surface. Our discussion is limited to the creep-like 
motion, when the driving velocity is small. Two regimes of motion are 
observed: stick-slip and steady sliding. The stability of the steady sliding 
depends on the model parameters. Numerical and analytical results show 
a transition between two regimes: the system passes directly from uniform 
to stick-slip motion. The calculations are performed also for two driven 
blocks. Then, the transition has the same character and it appears at the 
same point. 
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1 Introduction 

Friction between solids is a complex problem, first investigated in 1699 by Amon- 
tons In 1781, Coulomb formulated his laws, known as Amontons-Coulomb 
laws. Suppose that a solid block pressed to a flat surface by a normal force W , 
and the contact surface is of nominal area S. The static and dynamic friction 
coefficients /i^ and as the proportionality coefflcients between the friction 
forces Fs and Fd and the force W . Fs is the force necessary to move a standing 
block. Fd is the friction force during a uniform motion. The laws state that 
both friction coefflcients do not depend on W and S, and that fid < iJ-s- These 
laws are widely used until now 

In recent years, interest on friction has been revived because of its possible 
relevance for earthquakes and of its obvious relevance for technology 0. 
Detailed experimental studies revealed subtle effects on nanoscopic scale, and 
several phenomenological models have been formulated ^|Sl|Hl|7j. In particular, 
the creep regime is of interest, which can be considered as an intermediate stage 
between motion and rest. The idea is that during a very slow motion, the 
interface between solids preserves to some extent the information on history of 
contact points. For this regime, the model of Rice and Ruina "6 is considered 
to be appropriate [SJ |5] . Some variation of this model was formulated in T' . 

In these models, the physical system is a block connected to a driving mech- 
anism by a spring. The mechanism drives the block with a constant velocity. 
The uniform motion of the block is not a unique solution. An alternative is the 
so-called stick-slip motion, when the block velocity is varying periodically from 
zero to some maximum, then back to zero and so on. There is an experimental 
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evidence that once the uniform phase loses its stabihty, the motion becomes 
osciUatory, and the amphtude of the osciUations increases gradually from zero 
when the control parameter departures the transition ^[71 El E| • This sug- 
gests an appearance of the supercritical Hopf bifurcation ^Jj. However, as it 
was pointed out in |12| . further corrections to the Rice-Ruina equations are 
necessary to reproduce the gradual increase of the amplitude of the oscillations, 
present in the Hopf bifurcation. Without these corrections, the transition leads 
directly from the uniform motion to the stick-slip effect. In other words, the 
Rice-Ruina equations are not generic; a slight modification of the model leads to 
finite oscillations, and the Hopf bifurcation is observed. Examples of such mod- 
ifications are described in |12l llHj ; more general discussion of possible memory 
effects in the problem of friction can be found in ^3]. However, the Rice-Ruina 
equations are in most cases sufficient for a convenient description of the uniform 
motion in the creep regime. 

In realistic situations, the number of contact points between the moving 
surfaces is larger than one, and it depends on the normal force 00 El- This 
is an indication of a gap between theory and experiment on friction. It is 
clear that a more appropriate and general model should contain the number 
of contact points varying in time, with a distribution of elastic forces between 
them. However, results of such a model are expected to depend on numerous 
parameters, which vary not only from one sample to another but also in time. 
In particular, the number of contact points and their elastic constants are hard 
to be controlled. 

The aim of this paper is to compare the instability of the uniform motion for 
one and two blocks. In this approach, the blocks are equivalent to the contact 
points. Our main goal is to prove that the instability of the uniform motion oc- 
curs for the same values of the model parameters for one and two blocks. When 
discussing this result, we are faced with almost all above mentioned difficulties 
in the interpretation. In reality, the case of two contact points is probably as 
rare as the case of one. However, it seems to us that at least the direction is 
proper. We hope that our analytical and numerical results provide a basis for 
more extensive search. 

In subsequent section we show that for one block the analytical conditions for 
the Hopf bifurcation are not fulfilled. The calculations are supplemented 
with some numerical results, reported in Section III. In Section IV we prove 
analytically, that the transition point is the same for the case of one and two 
blocks. Again, the simulations support these results. Last section is devoted to 
final conclusions. 



2 The equations of Rice and Ruina 

The equations of motion can be formulated ^[^1 as follows 

§iVt-.)^^.o + Bln§-+Aln^^ (1) 

4>^l-f (2) 

where x — Vt is the block position with respect to the driving mechanism, V 
is the driving velocity, K is the spring constant, W is the normal force, /io is 
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a reference value of the friction coefficient for steady sliding at some velocity 
Vq. During sliding, the microcontacts are refreshed, on average, after a distance 
Dq. The state of these microcontacts is described by the variable 4>, which 
interpolates between the time of stick for the block sticked and Dq/V for steady 
sliding. Finally, <I>o = Dq/Vq and A, B are unitless material constants. We note 
that _B > A in the experimental data 

The equations can be transformed to an autonomous form. Denoting {x — 
Vt)/Do = a, $/$o = Vot/Do = r, V/Vo = to, exp(-^oM) = 7, KDo/W = 
K, we get unitless equations 

d = — w + 70~^ exp[— — a] (3) 

0= exp[--Ja] (4) 

where the time derivative is over r. These equations can be furter simplified by 
a change of variables: —l/ui exp{—Ka/A = x, (j)^^^^ — y, where b = B/A — 1. 
Introducing a parameter fx = b — k/A, and renormalizing time once more (r — s- 
7t) we get the equations in more algebraic form 

X = {b — ii){mx — x^y) (5) 



y^{\ + b){xy'' ^r^y^+ — ) (6) 

where m — lo/^ > Q, ri — I/7. At the fixed point, where the block velocity is 
constant and equal to the driving velocity, x — Q and y = 0. There, 

^ = v{-f (7) 
m 

. (1+&) 



y 

The determinant of Jacobian J at the fixed point is (& — ^)m? , and the trace 
is ^.m. This means that the stability of the fixed point is lost when the control 
parameter /i becomes positive ^T]. Near this point, Det{J) > 0. 

In due course, only first derivatives over fi will be calculated at /i = 0. Then, 
writing down the eigenvalues of J we can neglect terms proportional to fi^ . 



A = m(|±^y&^) (9) 
The transformation to the Jordan form leads to new variables ^,■0 

e = + (10) 

(11) 

V 

Now, the equations of motion are 

(1 + b) 2 b ?] b 

2b rj rj 
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Figure 1: Velocity of the single block versus time: (a) stick-slip (b)steady slip 
motion regime. 



^ = mVb{l - _ ^b^b)^/ - (1 + bM-)^^'+Th (13) 

2o 7] rj 

These equations can be written in short as ^ = /(CiV')i V' — ff(C)V')- The 
condition for the presence of the Hopf bifurcation (llj is 

+ /0V') ~ .9«v(5« + 9i>'ip) ~ fa9i£. + Ui>9i>i>] 7^ (14) 

However, direct calculations for Eqns. (12,13) at the fixed point and /i = lead 
to the result a = 0. 

As noted above, the phase of steady slip ceases to be stable when ^ = 
B/A — 1 — KDq/ (WA) becomes positive. This can mean in particular, that the 
spring constant K decreases. This is in accordance with the phase diagram in 
the creep regime, observed experimentally [Tl [71 lUl [TUl [T3] . 



3 Numerical calculations for one block 

The stability of the solution of the equations of Rice and Ruina is checked 
numerically at both sides of the transition. We applied the Runge-Kutta method 
of 4-th order. The result is shown in Fig. 1. The parameters of the calculation 
are: B = 0.08, A = 0.03, ^Io = 0.4, V = OA^im/s, cjio = 1.0, L»o = O.l/im, 
K/W = 0.054/.tTO"^ for the steady shp and K/W = 0.046/xm"^ for the stick- 
slip. The initial value of $ is slightly (0.01) different from value characteristic 
for uniform movement with velocity equal to driving mechanism velocity V . The 
initial values of block's position and velocity are equal to values characteristic for 
uniform movement. As we see, either the steady slip or the stick-slip is observed, 
without an intermediate phase with a continuous rise of the amplitude of the 
oscillations. Such a continuity is expected at the Hopf bifurcation. 
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4 The case of two blocks 

In this case we rewrite the equations (1,2) twice with new variables xi,x2,^i, $2 
and wc add a coupHng between the blocks by means of a new spring with 
constant ^2. Both driving springs have the same constants ki. The equations 
are 

/io + Bln Ji+^ln^ = ^{ki{vt-xi) + k2ix2-xi)) (15) 
Vq yy 

$2 Xo 1 

Ho + B\u— + A\u— = T^Aki{vt - X2) + k2{xi - X2)) (16) 

sPq ''0 

^1 = 1-^ (17) 

$2 = 1- -=— (18) 

Changes of variables, similar to the ones applied above, lead to autonomous 
dimensionless equations 

ai = -w + ^i ^ exp(-/io)exp(-^ai)exp(-^(a2 - ai)) (19) 

a2 = -w + 02 exp(-/;io) exp(-^a2) exp(-^(Q;i - ^2)) (20) 

= l-7(;ii"^exp(-/xo)exp(^i^ai)exp(^^(Q!2-Q!i)) (21) 

(^2 = 1- l(j)2~^ exp(-/io) exp(^i^Q!2) exp(^2^(ai - ^2)) (22) 
The eigenvalues of the Jacobian at the fixed point are 



Ai = -^[-Doki -AW + BW - ^/-AADokiW + {D^ki + AW - BWf] 



(23) 



A2 = [-Dofci -AW + BW + ^J-AADok^W + (£>ofci + AW - BWy] 



(24) 



^3 = TTJTrA-Doki - 2Dok2 -AW + BW 
2 AW 



^/{Doki + 2Dok2 + AW - BW^) - AAWiDoki + 2Dk2)] (25) 



^4 = TT^W [-^oA;i - 2Dok2 -AW + BW + 

2 AW 



y/{Doki + 2Dok2 + AW- BW^) - 4:AW{Doki + 2Dk2)] (26) 



5 



3 I 1 1 1 n — \ — I 0.1015 




Figure 2: Velocities of two blocks versus time: (a)slip-stick(b) steady slip motion 
regime. 

Two of them, Ai and A2, change sign at the same value of ki/W where the 
transition for one block occurred. At this point, A3 and A4 remain negative. 
This completes the proof that the transition point for two blocks coincides with 
the transition for one block. This is confirmed numerically, as shown in Fig. 
2, for the same parameters as above, and k2/W = 0.05/xm~^. In this case, 
the numerical plots reveal almost full synchronization of two blocks below and 
above the transition: the time dependences of the block velocities coincide. 

5 Conclusions 

We demonstrate that for two contact points, the uniform sliding ceases its sta- 
bility at the same value of the elastic driving force, that for one contact point. 
This means, that the elastic force between the contact points, represented here 
as k2, does not influence the stability of the uniform motion. A question arises, 
if this result could be valid for a larger number of blocks, i.e. of contact points. 
Such a generalization is of obvious interest; if true, it could release the limi- 
tation of the whole approach, which in this case could apply to real surfaces 
with multiple contact points. However, we can only state that the question still 
remains open. 

For two blocks, our numerical results provide a demonstration of a synchro- 
nization of the stick-slip motion. However, we know that such a synchronization 
can depend on the initial conditions, and therefore it cannot be treated as generic 
citemy. This is true in particular if one tries to generalize the results for the 
larger number of contact points. 
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